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Complex Path Integrals and the Space of Theories 

D. D. Ferrante, G. S. Guralnik, Z. Guralnik, C. Pehlevan 

The Feynman Path Integral is extended in order to capture all solutions of a quan- 
tum field theory. This is done via a choice of appropriate integration cycles, 
parametrized by M e SL(2, C), i.e., the space of allowed integration cycles is re- 
lated to certain Dp-branes and their properties, which can be further understood 
in terms of the "physical states" of another theory. We also look into representa- 
tions of the Feynman Path Integral in terms of a Mellin-Barnes transform, bringing 
t-h the singularity structure of the theory to the foreground. This implies that, as a 

sum over paths, we should consider more generic paths than just Brownian ones. 
Finally, we are able to study the Space of Theories through our examples in terms 
of their Quantum Phases and associated Stokes' Phenomena (wall-crossing). 
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1. Introduction 



This paper has the goal to extend our previous work, [1], bringing to the foreground some 
issues of Quantum Field Theory that have not quite found their place yet. In particular, it aims 
to introduce quantum phases as the analogous object to chamber walls (wall-crossing and 
Stokes' phenomena [ I]), in order to study the "Space of Theories" via suitable generalizations 
of the [Feynman] Path Integral. 

We will start by considering the quantization procedure as defined by the path integral and 
its respective algebra of observables, 



where is the Action of the theory with S representing a possible surface term needed in 
order to establish the self-adjointness of the Laplace-Beltrami operator in the kinetic term of 
the Action, interpolating among the possible boundary conditions suitable for the problem at 
hand, and t represents the possible coupling constants of the theory; t = t [if], i = 1, . . . , n 
are operators forming an algebra in an appropriate Hilbert space; and 3~ represents the Fourier 
Transform (with respect to the source J), and the path integral and the trace of the operators 
are cast as suitable Fourier transforms of the Feynman weight, e lS ^ [(/?] . 

In order for the Fourier transforms above to be well defined, the integrand and the kernel 
have to be single-valued. To that end, we can consider their Complexification and look for 
suitable integration cycles labeled by t and E, C^, i.e., Riemann sheets rendering our problem 
well defined. If 6^ is ultimately R, we recover the above definitions in the sense of Fourier 
transforms. However, for non-trivial cycles we end up with a more general integral transform, 
one that preserves the symplectic structure of the problem at hand. In fact, we can ask our- 
selves what is the most general integral transform that preserves the symplectic structure of a 
problem. And, as we will see in sections 2 and 3, this is accomplished by the Linear Canonical 
Transform (a linear map of the phase space preserving the symplectic form). 

In this sense, it is clear that t and S will determine the analytic structure of our quantiza- 
tion procedure: in fact, the integration cycle 6^ ultimately determines the fall-off properties 
of the fields in our theory, thus relating the decaying properties of the fields with the analyt- 
icity of the integral transform at hand, [16, in particular, the Paley-Wiener theorem (and its 
various generalizations), describing functions through their analytic extensions in the complex 
domain]. Actually, making this relation explicit will motivate us, in sections 4 and 5, to cast 
the Path Integral as a Mellin-Barnes transform, bringing the role of the analytic structure of 
our theory to the foreground. 

Thus, in a sense that will become clearer later, the central concept will be the analytical 
structure of our quantization procedure via the Feynman Path Integral, including its source 
term, e lJ ^: both, the integrand and the source term, have to live in the same Riemann sheet, 
otherwise we get an ill-defined theory. Therefore, realizing that the source term reflects the 
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analytic structure of the Feynman weight (after all, J and if are related by a Legendre trans- 
form), we can already infer that if — > — i d 3 is not going to express the most general case: 
this will only be true in the situation where 6^ defines a Fourier transform; otherwise, the 
Legendre-dual relation between J and if will typically be more elaborate. 

As a side remark, we can note that in order for the path integral quantization to be consis- 
tent with its respective Schwinger-Dyson equations, the field, if, has to be determined by its 
respective quantum equations of motion. Therefore, given that t and £ determine the prop- 
erties (monodromy Stokes' phenomena, etc) of the quantum equations of motion, i.e., of the 
Schwinger-Dyson equations, the integrand and e lJif have to live in the same Riemann sheet. 
This is a trivial observation if the quantum equations of motion are of 1 st degree. However, 
they are typically of higher degree, meaning that there are multiple solutions, which we can 
label through t and S, thinking that each solution is associated to a quantum equation of 
motion with suitable values of t and S. 

Thus, in order to bring these issues to the foreground, we will generalize the path integral 
(thought as an integral transform in a general sense, or as a Fourier transform in a more 
specific one) in two different ways: as a Linear Canonical Transform (sections 2 and 3), and as 
a Mellin-Barnes Transform (sections 4 and 5). These are particularly useful, since the former 
aims to be invariant by canonical transformations, while the latter is relevant to make explicit 
the analytic structure of the theory, e.g., the modern subject of scattering amplitudes in N = 4 
Super Yang-Mills. In this way, we will be able to associate certain theories formulated as Linear 
Canonical Transforms with their respective Mellin-Barnes transforms, making clear how one 
framework is related to the other, i.e., we will be able to relate certain orbits of the SL(2,C) 
parametrizing the Linear Canonical Transform (as seen below) with appropriate Mellin-Barnes 
transforms, connecting the symmetries described by these orbits to the singularity structure 
made explicit in the Mellin-Barnes transforms. 

This brings us to the question of what kinds of paths are allowed within the context of the 
path integral. After all, taking the analytic structure of the quantization process, of the path 
integral, seriously, how do we know that more exotic paths (woven across several allowed 
and distinct Riemann sheets) are not contributing, despite cancellations due to the highly 
oscillatory nature of the integrand? 

In quantum mechanics, the usual story says that all paths are allowed, although the equiv- 
alence between the integral formulation (path integral) and the differential one (Schrodinger 
equation) is only established at the level of Brownian paths, where we think of Schrodinger's 
equation as a diffusion equation, analytically continued. We are interested in this connection 
because of its possible implications in Quantum Field Theory, where we deal, instead, with 
Feynman's Path Integral and Schwinger-Dyson's equations as the two sides of the quantiza- 
tion procedure. We want to study the connection between the dimension of the solution space 
of the Schwinger-Dyson equations and the number of vacuum states present in a physical 
theory: in quantum mechanics, the dimension of the solution space of Schrodinger's equa- 
tion determines the number of ground states a certain model has — as is the case with all 
differential equations and the dimension of their respective solution spaces: the dimension 
of the solution space (vector space of analytic functions) of a n th order linear homogeneous 
differential equation with analytic coefficients is n. How does the equivalent problem work in 
Quantum Field Theory? Does Feynman's Path Integral include only Brownian paths or does 
it include more generic paths? It is to make sense of questions such as these that we want 
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to study the connections between the integral and differential formulations of Quantum Field 
Theories. 

The first extension of this result was done in [2] and can be achieved if one consider slightly 
more generic paths, called Levy flights. In this case, the path integral becomes a Fractional 
Fourier Transform. The gist of this approach is that Schrodinger's equation becomes fractional 
(i.e., rather than differential operators one uses fractional differential operators), while the 
path integral now includes paths of fractal dimension. More to the point, while Brownian 
motions have a scaling relation of the form x 2 oc t and fractal dimension df™™™ n = 2> Levy 
flights scale as x a oc t and have fractal dimension d f ^ al = a, with < a ^ 2. This means 
that, for example, in the case of a nonrelativistic particle, its dispersion relation is given by 
E = D a p a , where D a is the generalized fractional quantum diffusion coefficient. If a = 2, we 
have that D 2 = l/2m and E = p 2 /2m, which is the case for Brownian paths. 

The connection between this analysis in terms of scaling and dispersion relations, and our 
previous argument, is through the fall-off properties of the field content of our theory. Roughly 
speaking, the integration cycle 6^ ultimately determines these decaying properties (Stokes' 
phenomena, monodromy of the quantum equations of motion, etc). And, as such, 6^ is re- 
lated to the scaling and dispersion relations above, determining whether they are isotropic 
or anisotropic (and the respective scaling factors). Finally, the subject of anisotropic scaling 
relations has been studied before in Hofava-Lifshitz gravity, where the scaling goes as x oc t z , 
meaning that a ~ 1/z, and considering the Lorentz symmetry as emergent when the theory 
flows to some value of z ~ 1/a (e.g., z — 3 = 1/a), [6]. Unfortunately, the relation between 
anisotropic scaling, fractality and integration cycles of the path integral does not seem to have 
been analyzed before. 

However, this is not the most general case one can consider. In fact, this is just one among 
several possibilities. The bottom line is that the path integral can be parametrized by SL(2, C), 
where different orbits of the group correspond to distinct possible transformations: from a 
Fourier and Fractional Fourier, to a Linear Canonical Transform. 

Furthermore, we will use the expression of the path integral as a Linear Canonical Trans- 
form in order to introduce a representation of the path integral as a Mellin-Barnes transform, 
given by Fox's H-function. This is possible, as we will see later, because of a result stating that 
all stable probability densities can be represented in terms of Fox's H-function, [9]. The advan- 
tage of such a construction is that we make explicit the analytic structure of our theory: the 
monodromy of the integration cycles defining the H-function representing the path integral 
will determine the ramification of the theory at hand. The study of ramifications is generally 
described by Picard-Lefschetz theory. 

Thus, we tie in with our previous construction in different ways. Firstly, this approach has 
been used before, namely in [2], in order to express the fundamental solutions for a frac- 
tional diffusion process in terms of an appropriate H-function. Secondly, the properties of 
the SL(2, C) parametrization of the Linear Canonical Transform are reflected in the proper- 
ties of the H-function: given that we explicitly have the analytic structure of the theory, it is 
not difficult to see the results of certain analytic continuations. Finally, different H-functions, 
representing different path integrals, will be associated to different stable probability distribu- 
tions and, as such, with different allowed paths in the path integral, giving further evidence 
that we need to include more exotic paths than just Brownian ones. 
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At this point, a note is in order: the above results are not dependent on the particulars 
of our example (s). For instance, if could be R- or C -valued, it could also be Matrix- valued 
(or, more generally, Tensor-like in character), or even f£ ie algebra-valued (including graded 
ifie algebras, i.e., SUSY): the exact same results would be true (with minor modifications, 
e.g., tracing out the matrix degrees-of-freedom, or defining the Fourier Transform over the 
appropriate group, etc). 

The outline of this paper is as follows: in section 2 we will introduce the Linear Canonical 
Transform, and in section 3 we will describe the path integral as an LCT; in section 4 we will 
introduce Fox's H-function and some of its generalizations, while in section 5 we will describe 
the path integral as a Mellin-Barnes transform via the use of the H-function; in section 6 we 
will develop several examples supporting our previous constructions, and in section 7 we will 
suggest how this approach can be used to deal with higher-dimensional systems. Finally, in 
section 8 we will make a synopsis and hint at possible further developments. 



2. Linear Canonical Transform 



Canonical Transforms are those transformations which map the operators Q and P into linear 
combinations of themselves. The Linear Canonical Transformation is the linear map which 
preserves the symplectic form: it is understood as a family of integral transforms, parametrized 
by SL(2, E), generalizing several classical transforms, such as the Fourier, the fractional Fourier, 
the Laplace, the Gauss-Weierstrass, the Bargmann, and the Fresnel transforms. 

One of the Fourier Transform's main properties is that the operator Q (multiplication by 
the argument: (Q/)(q) = q/(q)) is mapped into the operator P (differentiation: (P/)(q) = 
— i ^p), and vice-versa with a minus sign. Such behavior is analogous of a rotation by n/2 
in the Q-P plane, the Phase Space. However, we will look for linear operators C which map Q 
and P into linear combinations of each other, 



where a,b,c,d GK for now, with the caveat that [Q 7 , P 7 ] = C [Q, P] C' 1 = i t, implying that 
ad — be = 1. 

We readily see that the Fourier Transform corresponds to a = = d, b = 1 = — c, and 
that the identity transformation corresponds to a = 1 = d, b = = c. So, we will label the 
transform operator as C M with the unimodular matrix 



q , = CQC~ 1 = dq-b¥; 
P / = CPC" 1 = -cQ + aP; 




detM = 1 . 
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The operator C M is linear and has the following properties: 

(C M /)(q)=/ M (q); 
(OmQ/Xq) = (CmQC^ C M f)(q) ; 

'd / M (q) 
dq 



= d(qf M (q))+ib^- 
(C M P/)(q) = (CmP^ 1 C M /)(q) ; 



= -c(qf ^)- ia \-^-j 



Integral Transform and Kernel To concretely realize the operator C M l et us use an integral 
transform with the kernel C m {q' ,q), 



f M (q') = (C M fXq / ) = 



f(q)C M (q / ,q)dq; (1) 



Combined with the properties previously derived, we have the following expression for the 
kernel C M , 

CM,q) = i^e^^ 2 " 2 "^^/" . (2) 
V2nb 

For Meal parameters, (2) oscillates strongly in the limits where \q\ — > oo and |q'| — > oo but with 
a fixed modulus at 1 / y/2n\b\. 
The inverse of the C M transform (1) can be shown to be given by 



/(q) = (C w 7 M )(q) = 



f M (q')CUq,q')dq'. (3) 



The composition of two such transforms, C Ml and C M , is given by the transform whose 
matrix parameter is the product of the constituent transforms, i.e., Cm 2 m 1 - I n Group Theory, 
the C M are a ray representation of SL(2,IR) (2x2 unimodular matrices), also known as the 
metaplectic representation (Weil). 

The identity C M transform corresponds to M = 1, as expected: lim M _^ 1 CmOz'jQ) — <5(q — qO- 
From this result, it is possible to show that all lower-triangular transformations M, where 
b = 0, have the form, 

Jcq 2 /2a 



(C M (b=o)/)(q) = — —f{q/a) 



e 

\fa 

These are called geometric transformations, once they consist of dilations by a and/or multipli- 
cation by an oscillating Gaussian. 
As for the inverse canonical transform, C M -i(q',q) = C^q,^), we just have to note that, 



M -1 = 



d e~ in b 
e lK c a 
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where we have to watch out for the appropriate sheet in the complex b -plane. 

Thus, we have showed that the set of LCTs C M forms a group of unitary transformations, 
SL(2, R); i.e., LCTs are a family of integral transforms parametrized by a 3-manifold that can be 
understood as the action of SL(2,IR) on Phase Space, the field-source plane. The name comes 
from the map preserving the symplectic (or multiplectic, for the case of multi-dimensional 
transforms) structure, canonical transformations, once SL(2,M) can be understood as the sym- 
plectic group Sp 2 . Thus, the LCTs are linear maps in the coordinate-momentum domain Q-P 
(i.e., phase space), preserving the symplectic form. 

The LCT, parameterized by M e SL(2, M), M = ( a c \ ) with a d - b c = 1, can also be defined 
more straightforwardly as follows, 

n oo 

JZi e inJ 2 d/b p ^ e in V 2 a/b e -2ni V J/b d(p . b ^ ; 

J —oo 

^e inj2cd F[Jd]; b = 0; 

where the prime denotes those elements where b = 0: M 7 = (^ °) G SL(2,R). 

Let us describe some of the special cases generalized by the LCT, i.e., physically relevant 
cases characterized by symmetries described by orbits of SL(2, C): 

Parabolic Subgroup or Dilation Transform Corresponds to a dilation by a factor A (and 
are a family of parabolic elements, represented by a = e _/3/2 , [3 e R), 

«-(;:)=(::-)■ 

Fresnel or Gauss-Weierstrass Transform The Fresnel Transform corresponds to shearing 
by A (cutoff), while the Gauss-Weierstrass "diffusive" transform is its analytical continu- 
ation, A —> A = —i b: 

Gaussian Transform Multiplication by a Gaussian of width —i w, 

»--(: !)-(^ 9- 

Elliptic Subgroup or Fractional Fourier Transform Corresponds to a rotation by an arbi- 
trary angle co; they are the elliptic elements of SL(2,M), represented by: 

fa b\ f cos co sin co 
M FrFT = A \ = \ • 

hrM yc dj y—smco cosco 

Fourier Transform Corresponds to a rotation by 90°, which is a special case of the Fractional 
Fourier Transform for M FrFT (co = tt/2), with the corresponding integral kernel given by 
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C M ^ = e- i7r/4 M FT , where 



Hyperbolic Subgroup 



a M ^ cosh co — sinhco 



^ HS I c d i I — sinh co cosh co 



Laplace Transform The bilateral Laplace Transform can be obtained once we extend into 
the Complex domain, i.e., SL(2, C), 



Bargmann Transform This is a particular case of the Complex LCTs, analogous to the Fourier 
Transform for Meal LCTs, where we use —i = e~ inl2 



For transformations with a^O (i.e., excluding FTs and special cases of the FrFTs), we can 
show that every C M £ SL(2,IR) can be written as the product of a GT, a DT and a GWT, as 
follows: 

When involving the Fourier Transform, we have that, 

As expected, the generalized uncertainty relation we get from and A^m, respectively the 
dispersion of / and f M , is given by: 

Af AfM ^ — . 

Finally, we notice that the transform kernel, C M , plays the role of the Dirac 5, in the sense 
of closing the relation of completeness. This is called the Coherent-State basis. 
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Hyperdifferential Operator Realization We want to find a differential operator in the 
following context: given a 1-parameter family of integral transforms, 



(C MM f)(q') = 



(4) 



including the identity for t = 0, i.e., M(0) = 1 and /(q 7 , 0) = / (q0> we want a differential 
operator H that can be written as, 



(C MW /)(q) = exp(i t ^0/(q) = J] L_J_ /(q) 



n=0 



n! 



Operator of this kind appear mainly in connection with the time evolution of the wave and 
diffusion equations, translations and dilatations. 
Formally, we can write 



CM(T)(q / J q)/(q)dq 



T=0 



The operator H generates the 1-parameter integral transform family (4). 
Of particular interest are the following 1-parameter subgroups of SL(2,R): 



Parabolic Subgroup or Dilatation Transforms = ^(QP + PQ), 



M DT (e- T ) = 



r T o 

e T 



dr 



1=0 



= l5(q-q')-q'^5(q-q'); 
= (\- q Ta) C ^ {q, ' q \ =0 



Fresnel or Gauss-Weierstrass Transforms 9f 



•GW _ 1 m>2 



M gwt (t) - 



— Cm gw1 "(<2 ,q) 

O T 



1 -T 

1 

1 d 2 



T = 



2 3q : 



C m g\nt (q',q) 



T=0 



Gaussian Transform = 



— 1 rT^2 



M gt (t) = 



1 

T 1 



i 



T=0 



= -q 2 C M GT(q / > q) 



T=0 
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Hyperbolic Subgroup ^C HS = - (P 2 - Q 2 ), 



M H5 (t) = 



T=0 



cosh t — sinh t 
— sinh t cosh t 

i A 5 2 

2 



T=0 



Elliptic Subgroup or Fractional Fourier Transform = - (P 2 + Q 2 ), 



IM (t) = 



cos t — sin T 



<7T 



T=0 



sin T COS T 

i f 5 2 

2 



+ q 2 GvMq^q) 



T = 



Again, if we allow for Complex parameters in the expressions above, we obtain the group 
SL(2,C) of Complex 2-dimensional special (unimodular) linear transformations. 

We now see that linear canonical transformations are generated by all operators constructed 
out of quadratic expressions in Q and P. Further, we have the following objects, 

1. Integral Transforms C M ; 

2. Hyperdifferential Operators exp(i t?T); & 

3. 2 x 2 matrices M. 

For every element in one there are corresponding elements in the other two, and this corre- 
spondence is preserved under composition, sum and multiplication. 



3. Path Integral as a Linear Canonical Transform 

One of the more common quantization techniques is via the Feynman Path Integral formu- 
lation, generalizing the Action Principle. Furthermore, the Path Integral allows us to easily 
change coordinates between different canonically conjugated descriptions of the same quan- 
tum system. However, it is rather unfortunate that most treatments of the Path Integral do 
not do it full justice, e.g., a priori assuming that its measure is M-valued, and thus missing 
on possible extra solutions to the quantum equations of motion and extensions of the theory 
(c.f., [7, 8, 10, 16], roughly summarized in appendix A). 

Historically, defining the Path Integral's measure, in any way (ranging from more mathe- 
matically rigorous approaches all the way to the more formal and heuristic ones), has always 
been the Achilles' heel of the Path Integral quantization scheme. What we will show below is 
a generalization that brings to the foreground the Legendre-dual relation between the field, ip, 
and its source, J. Once again, note that <p can be quite general: M- or C-valued, Vector- or 
Matrix- or Tensor-valued, and even ifie- or graded ifie Algebra-valued (SUSY): the neces- 
sary modifications, e.g., tracing matrix degrees-of-freedom, do not affect our discussion: there 
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is no loss of generality in the results below dealing with scalar fields (spin bosonic fields) or 
DO-branes. 

The idea of generalizing the Path Integral not only as an integral transform, but as a general 
Linear Canonical transform, explicitly bringing forth the symplectic structure of the theory is 
motivated by several reasons. In what follows, we can consider a "finite" (discrete) version 
of the Path Integral in several different but analogous ways: we can think of our theory as 
having an IR- and an UV-cutoff (including defining the theory via an appropriate Vertex Op- 
erator Algebra that implements the desired IR- and UV-cutoffs), or we can consider different 
discretization schemes (discrete differential forms, finite exterior calculus, Matrix models — 
BFSS/IKKT, infinite momentum frame: N x N matrices describing N Dirichlet particles — , 
naive lattice discretizations, etc). In these cases the Path Integral as well as the questions and 
motivations below are well-defined and tractable. 

Nonlinear Stokes' Phenomena & Path Integral Convergence: We want the Path Inte- 
gral to be well-defined in a certain sense, even if it is mathematically loose (and only sat- 
isfies the physicist's level of rigour). One of the ways to achieve this is through studying 
the asymptotics of the oscillatory Riemann-Hilbert problem associated to the nonlinear 
Schwinger-Dyson equation describing the system at hand: this means that the fall-off 
behavior of the fields of our theory will ultimately determine their algebra and quantiza- 
tion procedure. One of the methods used to attack such problems is the steepest descents 
(also called 'saddle-points' or 'stationary phase'), which identifies the contribution from 
the dominant saddle-point: once the integration contour is appropriately deformed, the 
original integral is extended (by a change of variables) to the entire Meal y-axis. In 
this sense, we will consider the Path Integral over a Complex integration cycle (a multi- 
dimensional integration contour), that we deform in order to obtain a Path Integral over 
the Meal if -axis (which is the usual way they are presented to us). Moreover, we want 
to associate to each integration cycle a Matrix label/parameter whose job is to lift' the 
Feynman Path Integral to an integral transform called Linear Canonical Transform, [3]. 

A previous generalization of the Path Integral along these lines has already been done, 
albeit only using the Elliptic subgroup of SL(2,M): it re-interprets the Path Integral as a 
Fractional Fourier Transform, [2]. Unfortunately, this is but a subset of all that can be 
done, although we can already see some of the features of the general case: the Elliptic 
subgroup of SL(2, M) consists of rotations parametrized by an arbitrary angle, 



which means that we are performing a rotation of co in Phase Space (if -J plane). This 
means that the Fractional Fourier Transform, i.e., the Path Integral, has real and imagi- 
nary parts, and its imaginary part only vanishes for special values of co. 

However, this rotation is limited by the integration cycle defining the Fractional Fourier 
Transform (viz. Path Integral): we cannot cross the poles of the integrand being trans- 
formed, given that the integration cycle delimits the allowed regions. The boundary 
between two allowed regions is subjected to Stokes' phenomena, also known as wall- 
crossing phenomena in this setting, [12]. 




— sm co cos co 
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These Stokes' lines are curves in the Complex ip -plane, where the integrand of the Path 
Integral decays without oscillations; analogously, anti-Stokes' lines are curves where the 
integrand oscillates without change in amplitude. Upon crossing the Stokes' lines, the 
integrand jumps discontinuously; equivalently the integrand becomes more and more 
wildly oscillatory as it approaches an anti-Stokes' line. The steepest descent path is 
relevant for topological reasons (more about this on the next item below), i.e., we need 
to know which saddle-point lies on the path. Lastly, we should remind ourselves between 
the connection of Stokes' phenomena and the study of hyperasymptotics: the analysis 
of divergent asymptotic series, where key parts of the answer lie 'beyond all orders', 
because the desired features are exponentially small in the perturbation parameter, [8]. 
This should gives us a clear connection between the perturbative series for our model, 
Stokes' phenomena, and wall-crossing. 

Along these lines, we can define each "branch" (wall chamber) of the theory determined 
by S T [ip], as bounded by appropriate Stokes' lines (codimension one manifold). Thus, 
we can translate the quantities computed in one "branch" (wall chamber) into another 
using the discontinuous jump determined by crossing Stokes' lines: this enables the com- 
putation of the discontinuous change of quantities such as integer invariants, space of 
BPS states, etc, across walls of marginal stability. In this sense, maybe we can borrow 
the definition of quantum phases, which originally means to indicate different quantum 
states of matter at zero temperature: the idea is that varying parameters of the theory, 
quantum fluctuations can drive phase transitions. In our current case, the integration 
cycle and associated Stokes' phenomena essentially determine the allowed range of con- 
tinuous values of the coupling constants. It is at the point where the coupling constants 
suffer a discontinuous leap (along a Stokes' line) that a quantum phase transition hap- 
pens. Thus, a given theory S T [(p] may have more than one quantum phase associated 
to it, provided we find more than one integration cycle rendering the Path Integral well- 
defined. 

Here we depict a "single theory", 
S T [(/?], that has three quantum phases 
for different values of the couplings: 
T — T i> T 2> T 3- Separating each set 
of allowed couplings are Stokes' lines 
(codimension one manifolds). In this 
sense, we can say that this theory has 
three quantum phases (wall chambers), 
such that crossing a Stokes' line (wall- 
crossing) is analogous to a quantum 
phase transition. 

Different Quantum Phases are associated with distinct integration cycles that are bounded 
by the appropriate Stokes' lines. This yields distinct values for the matrix parameter of 
the Linear Canonical Transform, M e SL(2, C), as shown in the examples below. Further, 
because each quantum phase is delimited by Stokes' lines, crossing from one quantum 
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phase to another implies the picking up of a Stokes phase factor, i.e., a finite contribution 
that accounts for the discontinuity between the two quantum phases, and is related to 
the monodromy of the Schwinger-Dyson equations. 

Finally given the nature of the Linear Canonical Transform as generalizing several dif- 
ferent kinds of integral transforms, it also represents an extension in the allowed classes 
of paths present in the Path Integral (viewed as a sum over paths) : on top of Brownian 
paths, one should also consider more generic cases, such as Levy flights, and so on, [2]. 

Symplectic Structure & Complex Path Integral: We can consider the above as the prob- 
lem of quantizing a symplectic manifold, namely the Phase Space of the particular the- 
ory we have at hands: we choose a Lagrangian splitting (a maximal isotropic subspace, 
a choice of coordinates and momenta, i.e., Phase Space P ~ L ® 1/ where L and 1/ are 
isotropic subspaces and dim R L = dim R V = V2 dim R P: this is tantamount to choosing a 
polarization for the phase space at hand), and the integration cycle for the Path Integral 
is the extra information provided towards integrability. 

In this sense, this is very reminiscent of the Bargmann-Segal representation (cf. Bargmann- 
Fock or holomorphic representation) equipped with a Kahler polarization and compati- 
ble affine and Complex structures, [17]. This is reasonable given that the Linear Canon- 
ical Transform also generalizes the Bargmann transform, as we saw in the previous sec- 
tion, with the matrix parameter being M BT = ^ ( } t \ l ) and —i = e~ inl2 . 

The relation between Phase Space quantization (also known as Weyl or deformation 
quantization: focus on the multiplicative structure of quantum observables), geometric 
quantization (focus on the representation of observables), and the A-model in string 
theory is discussed in [17, 18]. Essentially, the basic idea is to start with a symplectic 
manifold and choose a suitable complexification for it. There are three points of view to 
this subject: 

1. Complexify the phase space (i.e., the original symplectic manifold) to a complex 
symplectic manifold that has a good A-model with respect to the symplectic struc- 
ture given by the imaginary part of the complexified one; 

2. Start with a complex symplectic manifold and pick a suitable coisotropic brane, 
assuming it to be an A-brane with respect to the symplectic structure given by the 
imaginary part of the original complex one; & 

3. This is the more natural one from the viewpoint of 2-dimensional TFTs, bringing to 
the foreground the A-model of the complexified phase space: there may be many 
inequivalent choices of A-branes, corresponding to different choices of a complex 
structure. Thus, the same A-model can lead to quantization of the phase space in 
different symplectic structures. 

The point is to choose a suitable behavior at infinity for the functions that we can quan- 
tize. Generally, the condition chosen is that only functions of polynomial growth are 
allowed. However, the integration cycle that we choose for our Path Integral will ul- 
timately determine this fall-off behavior, [16]: understanding the Path Integral as a 
Fourier Transform, the idea is to relate the decay properties at infinity of its integrand 
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with the analiticity properties of the Path Integral (in a loose analogy with the Paley- 
Wiener theorem, [16]). This very same observation can be made through the study of 
the [global] boundary conditions of the Schwinger-Dyson operator (viz. Bethe-Salpeter 
equation, Ward-Takahashi identities, and anomalies). Thinking in terms of differential 
equations (local systems and Higgs bundles, [11]) makes the discussion about bound- 
ary conditions more straightforward, where different boundary conditions are related to 
distinct extensions of the respective differential operators, [10] (c.f., appendix A). 



In the general case, we have a G- 
bundle n : P — > M with fiber G. The 
different boundary conditions for the 
differential equations (local systems) 
are given by E : and S 2 . The fields of 
the theory are sections of the bundle. 
The red dots represent the points where 
the boundary conditions are satisfied. 
The Schwinger-Dyson operators (V) 
restricted to the global boundary condi- 
tions are V| E and V| E2 : these are two 
possible extensions of V, [10]. 

Historically, this study can be traced back to some early ideas by V I. Arnold, regarding 
Hermitian operators and its ground states, as well as the role played by solitons/instant- 
ons/monopoles (loosely referred to as "quasimodes"), [25]. Later, these were extended 
to symmetric operators with multiple ground states (c.f., appendix A, and [10]), and 
also to Complex structures on moduli spaces, [25]. 

Finally, the above can also be understood in terms of the Symplectic Field Theory de- 
veloped by Y. Eliashberg, [26], including its connection to Floer Homology. (Seiberg- 
Witten-Floer homology yields knot invariants formally similar to combinatorially defined 
Khovanov homology: these are closely related to Donaldson and Seiberg invariants of 
4-manifolds, as well as the Gromov invariant of symplectic 4-manifolds — solutions to 
the relevant differential equations, that we take to be the Schwinger-Dyson equations.) 

Effective Action & Schwinger's Action Principle: Expanding the motivations for this pro- 
posed generalization of the Path Integral, we now focus on the construction of the Leg- 
endre conjugation between <p and J, and explore its implications: 

if yr[j] = -i log 3T[J] ; and T[(p] = W[J] - J (p ; 

then ip^-ijjW[J]; and J^-i-^T[(p]; 

yielding r[<p] = S[<p] + \ logdet(^S[>]) . 



Typical fiber (G) 




» Base manifold (M) 



.^Boundary Condition 
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As long as the relation established by the Legendre transform between J and if is at 
most quadratic (in the expression defining the Effective Action, F, above), the Linear 
Canonical Transform is the most general integral transform preserving the symplectic 
(resp. metaplectic) form at hand, i.e., preserving the structure defined by the Legendre 
conjugation of <p and J. This implies that, in fact, we can even extend this construction 
from 1PI to 2PI, [15], and the Canonical Linear Transform can again be used to construct 
the Path Integral of the system, given that it still preserves the symplectic structure 
of the theory (once the 2PI construction saturates the quadratic relation between if 
and J, [15]). Speaking in terms of Deformation Quantization, the Linear Canonical 
Transform preserves the structure of the theory up to 0(R 3 ) deformations. 

Thus, considering the symplectic structure of the plane if-J, i.e., the theory's Phase 
Space, we can think of the case where (/?,JeCu{oo}, i.e., the field and source are valued 
in the Riemann Sphere: in such a scenario the symplectic structure is then preserved by 
the Mobius group, SL(2, C), which is the automorphism group of the Riemann Sphere. 
In the special case where the coefficients of a Mobius transformation are integers, we 
obtain the Modular Group, SL(2,Z), which is relevant for S- and T-dualities. 

Finally, with a bit of hindsight from (5) and (6), we can discuss the Legendre Transform 
construction of the Effective Action based on the Linear Canonical Transform expression 
for the Path Integral, in analogy to what is usually done when we treat the Path Integral 
as a Fourier Transform. Let us assume that b = — 2 n in the Linear Canonical Transform 
(cf. (5)), for if b = we only have trivial [free] dynamics (cf. (6)). Now, we have to 
discern between two cases: d = and d ^ 0. If d vanishes, we have that if — > —i5 Jy 
meaning that (if n ) = 5j2? M /(—2ni) n , implying (if) = 0. However, if d does not vanish, 

we get if — > j^jj an d (^) / 0- This gives us a direct way to account for symmetry 
breaking straight from the parametrization of the Path Integral given by M = (jj) G 
SL(2,C). 

In this sense, thinking of the Path Integral as a "sum over paths", we see once again 
that we are required to include more exotic paths than just Brownian ones, such as Levy 
flights, and so on: a generic Linear Canonical Transform, parametrized by a generic M e 
SL(2,C), will be characterized by different moments, (if n ), depending on the physical 
system at hand, ultimately describing distinct [stable] probability densities, [ ], that are 
associated to their respective paths and diffusion processes. Further, in anticipation of 
Section 5, it is worth noting that all stable probability distributions can be associated 
with a certain Fox H-function, [9]. 

Lastly, we can use the connection between stochastic processes and orthogonal polyno- 
mials in order to try and establish a Wiener-Askey polynomial chaos expansion of the 
Path Integral, [14], loosely along the following lines, 



i=0 

where {Z t } is the complete orthogonal polynomial basis of the Wiener-Askey chaos ex- 
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pansion, and £ is a random variable chosen according to the type of the random distribu- 
tion at hand, e.g., free fields, i.e., Gaussian random variables, yield Hermite-chaos, [14, 
cf. table 4.1]. 

The idea is to use appropriate Fox H-functions for the Z t % thus developing an "H-chaos 
expansion" for the Path Integral. This should become more clear in Section 5, when we 
venture into writing the Path Integral in terms of Fox's H-functions. 

Penrose Transform: Finishing the list of motivations for our generalization of the Path Inte- 
gral, we will borrow an analogy with the Penrose and the Penrose-Ward transform. The 
insight is that we are also trying to define some sort of Complex analogue of the Radon 
transform, where the straight lines we are looking for are not the light-cone as done in 
the Penrose Transform (massless fields, as in twistor theory), but the lines delimiting 
the Stokes wedge where a certain integration cycle is defined, i.e., the Stokes' lines of 
the integrand. In this sense, we would like to establish a parametrization of the space 
of integration cycles by the Matrix index, M £ SL(2, C), labeling the Linear Canoni- 
cal Transform, i.e., we would like to parametrize a certain space of compact Complex 
manifolds (the integration cycles) by a particular 3-dimensional submanifold of SL(2, C) 
determined by the allowed values of M. 

For massless fields, we already know that the Penrose Transform establishes a correspon- 
dence between two spaces, where one parametrizes certain compact complex submani- 
folds (cycles) of the other. Further, its non-linear extension, the Penrose-Ward transform, 
relates certain holomorphic vector bundles (on CP 3 ) with solutions of the self-dual Yang- 
Mills equations (on S 4 ). 

The idea is to construct a map between integration cycles (compact Complex manifolds) 
and solutions to the Schwinger-Dyson equations, given by the Linear Canonical Trans- 
form evaluated at a particular M e SL(2, C). 

Therefore, given all of the above, we will define our extended Feynman Path Integral to be 
given by the expressions below: 

3T M [j] = \[^ie inj2&lh jo[ip]e inLp2a/b e- 2ni ^ J/bf Dip ; b^O; (5) 

Sm'D/] = \fce inj2cd 0[Jd] ; b = ; (6) 

where M = { a c \) e SL(2,C) parametrizes our Path Integral, the prime denotes those elements 
where b = 0, i.e., M' = ( a c ° d ) e SL(2,C), and 0[>] is a certain functional of the fields. 

4. Fox's H-, and Generalized H-function Transforms 

The history of Mellin-Barnes transforms is intimately related to that of differential and dif- 
ference equations, as well as the development of hypergeometric functions: the robust un- 
derstanding of the asymptotic properties of functions, hyperasymptotics, Stokes' phenomena, 
monodromy and so on, [8]. 
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Fox's H-function was introduced to better the asymptotic control and to be the most gener- 
alized symmetrical Fourier kernel (a generalized Fredholm operator, [24]) involving Mellin- 
Barnes integrals. The H-function includes and generalizes several elementary higher, and 
special functions, ranging from the exponential, Airy Bessel, special polynomials, hypergeo- 
metric, Riemann Zeta function, up to Meijer's G-function, multi-index Mittag-Leffler function, 
polylogarithms, etc. This means that we can give a Mellin-Barnes representation to a wide 
range of functions. 

In particular, the H-function has recently found applications in a large variety of situations, 
ranging from diffusion all the way up to fractional differential and integral equations, and so 
on. It is worth noting that the H-function can be extended to have Matrix arguments (R- or 
C-valued), to include polylogarithms, and so on, [4,5]. 

The notation we will use for the H-function is as follows, 



H£ n 0O = H„T [« P A ; b 9 ,B Q ; z] = H„T 



PA 



PA 



(b q ,B q ) 



. Tjm,n 
H PA 



(a 1 ,A 1 ),...,(a p ,A p ) 
(b 1 ,B 1 ),...,(b q ,B q ) 



2ni 



c n ;= M+ i rd - b, - Bj 5 ) n; =n+1 n aj +a, S ) 



ds ; 



(7) 



where i = V— 1 , z ^ 0, z~ s = exp(— s(log|z| + i argz)), and argz is not necessarily the 
principal value, and an empty product is always taken as unity. There are three different 
integration cycles L, schematically depicted below, where red ticks are singularities of r(l — 
a i — Ajs) and green ones are singularities of T(bj +Bjs): 
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C r : goes from r — i oo to r + i oo, r e R, so 
that all singularities of T(b ; +Bjs) lie 
to the right of L l9 while all 
singularities of T(l — a ; — Ajs) lie to 
the left of it. 

C 2 : loop beginning and ending at +00 and 
encircling all the singularities of 
r(b ; +Bjs) once in the clockwise 
direction, but none of the singularities 
of r(l -CLj -AjS). 

C 3 : loop beginning and ending at —00 and 
encircling all the singularities of 
T(l — dj — Ajs) once in the anti- 
-clockwise direction, but none of the 
singularities of T(b ; +Bjs). 



The H-function is, in general, multi-valued, due to the presence of the z s factor in the 
integrand. However, it is one-valued on the Riemann surface of log(^). 
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For example, here are some functions expressed in terms of the H-function, i.e., the Mellin- 
Barnes representations of some functions: 



H, 



1,0 

0,1 



(b,B) 



= B- l z h ' B exp(-z 1/B ); 



(l-v,l) 
(0,1) 



= T(v) (1 + z y v = r(v) Mv; ; -z) ; 



(i-r,i) 

(0,1), (1-/3, a) 



= r( r )£^( z ); 



p,<z 



\a p ,A p = 1) 
L (b q ,B q = l) 



-.1,5 + 1 

J s+l,s+l 



0, 1 — a, . ..,1 — a 
0,-a,...,-a 



where r(z) is the Gamma function, p F q (a p ; b q ; 2) is the generalized hypergeometric function, 
E a p(z) is the generalized Mittag-Leffler function, G™^(a p ; b q ; z) is Meijer's G-function, and 
£(s, a) is Hurwitz's zeta function (that yields Riemann's zeta function as £(s, 1)). 

The H-function can be further generalized along similar lines as the construction done 
above, in order to include some functions which are not particular cases of Fox's H-function, 
e.g., Polylogarithms, including Polylogs of complex order, [5, Appendix A.5], and some further 
examples from Physics, [4] . 

For Polylogarithms and generalizations of the Riemann-zeta function, among others, we 
need to use one certain generalization of the H-function, called the H-function, H™'*(z). The 
H-function is defined as, 



H" 



(a ; -, GLj,Aj) l n , (a ; -, CLj) n +i tP 



2ni 



M L( b j^j)l,m,(fcj^j^j)m+l,q 

where #(s) is given by, 

n; =1 + p ]S ) uu rA ^ - a i - a ; s ) 



X(s)z s ds 



YY j=m+1 r^(i - bj - p jS ) n; =n+1 r(a, + a,s) ' 



(8) 



(9) 



where the contour £ iroo starts at r — i oo and goes to r + i oo, with r e R, separating all the 
singularities of T(bj — f}j s) from those of T(l — a ; - + a ; - s). For non-integer A ; - or Bj the poles of 
the Gamma functions in (9) become branch points, whose branch cuts can be chosen so that 
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the integration cycle can be distorted for each of the three £'s shown above as long as there is 
no coincidence of poles from any pair(s) of T's. 
Then, the polylog of Complex order v, L v (z), is represented as follows: 



L v (z) = Hl' 2 



2,2 



(0,U),(U,v) 
(0,l),(0,l,v-l) 



and a generalization of Hurtwitz's zeta function is given by, 



(10) 



<f>(z, q, rj)=H. 



- tjl.2 

2,2 



'(0,1,1), (l-^l.q) 
(0,1),(-T],l,q) 



— z 



k=0 



(rj + fc) q 



(11) 



The reason for our studying of the H-function is because it shows up when one generalizes 
from derivative operators to Fractional Differentiation, [2,5,9]. By extending the Path Integral 
to be represented by a Linear Canonical Transform, we have brought into play not only the 
elliptic subgroup of SL(2, C) — related to fractional differential operators — , but also the 
parabolic, the Fresnel and Gauss-Weierstrass, the hyperbolic, the Bargmann, etc, subgroups of 
SL(2,C), as already explained in previous sections. Moreover, as mentioned in the beginning 
of this section, the H-function is intimately related to the asymptotic behavior of the system it 
is describing, as well as its Stokes' phenomena and monodromy properties. In this fashion, if 
we can express the Path Integral as an H-function, we can condense all this information in the 
analytic structure of the H-function. This is the subject of the next section. 



5. Path Integral as a Fox H-, or Generalized H-function 
Transform 

Throughout the paper, we will also try and express the calculated Path Integral in terms of 
the H-function and its generalizations, [4]. The reason for this stems from the fact that the 
H-function is a Mellin-Barnes transform, which shows up in generalizations from derivative 
operators to Fractional Differentiation, etc, [2,5,9]. 

As such, we look for a representation of our theory in terms of a Mellin-Barnes transform, be- 
cause of its relation to asymptotic expansions (i.e., Stokes Phenomena and wall-crossing), and 
its closeness to Fourier Transforms: in this sense we would like to capture the properties of the 
theory straight from its singularities. We will consider the Path Integral as a hypergeometric 
integral transform, or simply an H-transform, [2, 19]. This view in terms of an H-transform is 
of particular significance given Fox's H-function's definition in terms of a Mellin-Barnes trans- 
form, that clearly brings the singularity structure of the particular H-function being used to the 
foreground. This shows that the Path Integral considered as an H-transform encodes all of the 
interesting and relevant singularity structure, which can be as simple as that of Polylogarithms, 
or far more intricate as we will show in the examples below. 

In the particular case of Polylogs, represented by the H-function, their recent importance 
for the computation of scattering amplitudes (maximally helicity violating amplitudes, MHy 
and Wilson loops in N = 4 SYM), their relation with Grassmannians and Cluster Varieties (viz. 
Stokes Phenomena), their relation to the Bloch-Wigner-Ramakrishnan-Zagier functions and 
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volumes of polytopes in AdS and mixed Tate motives, makes them a very relevant case to be 
studied. In this way, by expressing the Path Integral in terms of a generic H-function, we can 
more easily see what is induced by this structure. 

As an example of how the singularity structure of the H-function determines a theory we 
can look at [2], where the fundamental solution (Green's function) for a fractional diffu- 
sion process is given by a generalization of the Airy function, expressed in terms of Fox's 
H-function: H^[(l — 13/2,13/2); (0, 1); ip], where the limit /3 — > 1 recovers the Gaussian den- 
sity H};J[(l/2, 1/2); (0, 1); ip] = exp(-(/? 2 /4)/2 and the limit /3 2/3 recovers the Airy 
function, Hj;J[(l/3,2/3); (0, 1); ip] = 3 2/3 Ai(y/3 1/3 ). In fact, a stronger result holds: all 
stable probability densities can be represented in terms of Fox H-functions [9]. (A random 
variable is called stable (and has a stable distribution) if a linear combination of two indepen- 
dent copies of itself has the same distribution, up to translations and scalings. Examples of 
stable distributions include the Gaussian, Cauchy and Levy distributions.) 



6. DO-brane Examples 

We perform the analysis suggested above in three different examples, progressively showing 
more and more relevant details and properties of our extension for the Path Integral. 

The models we study include the free field, the cubic, and quartic interactions, revealing 
their different structure dependent on M e SL(2, C), i.e., using a 3-dimensional manifold as a 
parameter space. 

We would like to re-interpret our construction along the lines of a BFSS/IKKT matrix-model 
and consider its large-N limit, which is conjectured to be equivalent to M-theory. As such, 
rather than considering the problem from a bottom-up, constructionist, viewpoint, we would 
like to think about it using a top-down strategy: rather than solving a DO-brane model and 
considering infinitely many copies of it in order to build an extended model, we can think 
in "reversed" terms, where certain extended structures of the theory determine its properties. 
This has been named Brane Quantization, [18]. Our plan is to use our DO-brane examples 
in this section, considering the implications of the large-N limit of these models to certain 
structures in each one of them, e.g., the convergence of the Path Integral defining the model 
or, analogously, the boundary conditions of the Schwinger-Dyson equation each model obeys. 
For example, let us consider a Matrix-valued DO-brane model with a generic polynomial for 
an Action. Its Feynman Path Integral (FPI) and Schwinger-Dyson Equation (SDE) are given 

by, 

3f T [J] = Je itr ^=° A ^'V Jv D</> <oo; (12) 

e 

n 

tr (J] k A k (-i d,f~ l - j) ^ [J] = ; (13) 

k=i 

where if e C nxn is our Matrix-valued field and J e C nxn is its source, "Dip is the canonical 
integration measure in C nxn , C is an integration cycle rendering (12) finite, and t is a label to 
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remind us of the connection between 6 and {A fc }£ =0 (this will be further explored later), i.e., 
different sets of coupling constants will yield Actions that will be convergent under distinct 
cycles C, and t's role is to perform this bookkeeping. 

We could as well have written the above for a graded ifie Algebra-valued field, i.e., for a 
superfield 4> with superpotential W[$], 

M ab y abc 

ti, 5] = J e (L ° * a * c+ c - c0 e ir *« * a D<I>D4> < oo ; (14) 

e 

(M a b d$ - i y a bc d# d Z c - 5 a ) ^ T [Z] = ; (15) 
[M ab d ib -iy abc fy % - 5 G ) 3f T [5] = ; (16) 

where 3 a = 3 a + L a , i.e., the source is shifted by the linear term (ditto for its complex conjugate 
term). This motivates our treatment of the cubic Action on Section 6.2, i.e., cubic potentials 
are relevant to SUSY models, as well as to [non-Abelian] Chern-Simons theory. 

At this point we have some observations to make regarding the coupling constants of the 
model at hand, that analogously to the fields themselves, can be valued in M or C, or be Vector 
or Matrix valued, or ifie algebra or graded f£ ie algebra valued. 

The couplings can be understood as background fields, meaning that the renormalized, 
effective Action is constrained by symmetries (superselection rules), by its holomorphic prop- 
erties (Lee-Yang and Fisher zeros), and by its different properties in different limits (weak and 
strong coupling limits). Moreover, a holomorphic function (i.e., a section in an appropriate 
bundle) is determined by its asymptotic behavior and singularities, connecting us to he for- 
malism developed earlier, where the Path Integral is represented in terms of the appropriate 
H-function encoding its asymptotic behavior and singularities. 

Furthermore, we can implement desired boundary conditions for the Schwinger-Dyson 
equations in terms of a VEV(s) for the field(s) (at the boundary), (ip)\g M = const (flux quanti- 
zation), or via explicit terms in the action. See Appendix A and [10, 18]. 

So, considering the couplings as parametrizing a family of functions between two Morse 
functions (namely, the Hamiltonians described by the initial and final values of the couplings in 
a certain range), its degeneracies involve a birth (branching out) or death (merging together) 
transition of critical points. This study can be done via Cerf theory, which tackles stratified 
spaces labelled by the co-dimension of the strata, [25]. As such, in the examples below we 
will show plots to ilustrate these transitions at particular values of the couplings. 



6.1. Free Field 

In order to get our feet wet and highlight some of the properties that will be relevant later on, 
let us start with the simplest possible example: the free field. 

This DO-brane model has an Action given by S[<p] = /x y 2 /2, yielding the following FPI and 
SDE: 
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2f[J]=JV^ 2/ V vJ 2V<oo; 

e 

(-i fxdj — J) 2f[J] = => 2f [J] = a e ij2/2 " ; 



(17) 
(18) 



where a is an arbitrary constant, and 6 is an integration cycle. 

In order to cast it as an LCT, we will follow (5) and choose 0[i/?] = 1 and M = ( ±. % ). 

2 71 

This choice of M e SL(2,C) gives us the rational function determining the analytic structure 
of this transformation, 



Mt» = /„(</>) = 



-2 n (jjl <p + 2 it) 



(19) 



These are plots of f^ip) for different values of ju: 



ICii»+L.#l Lh-nh* F*#*m*nn Efelw* 






^'m = ' f! - -2 n (1 1+ ? ni'f 
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fiH--i.pl wiIm Ri«n«in Sphtr* 

Hu-'Lf] ■ -2 n ^ t 4 2 nj.i f 




Plots of f^(if) = ~ 2n ^^~ 4n t for /x = l,—l,i,—i # respectively. The plots on the right 
have BtCZ/^i^ii) as the x-axis and 3C//i=±i,±i) as the y-axis, with variations of color 
(phase 9) and brightness (magnitude p) to indicate z = p e l ° = ^=±1^; plots on the 
left are simply their wrapping over the Riemann Sphere. 



It is not difficult to see that these graphs are related to each other by a Mobius transform of 
the Complex <p -plane. 

For the special case where /x = i we recover a well known identity about Fourier Transforms: 
the Hermite polynomials are its eigenf unctions. This can be seen in (17) simply by noting that 
the th order Hermite polynomial is given by He (y) = 1, as defined in (31) and (32), and 
&[e~^ /2 He n (y)] = {i) n e~ J /2 He n (J), taking n = for this special case of ours. This result 
can be extended to other values of with a change of variables given by, if — > y/Tjl if and 

J — > J I yfiji, provided we stay within the same Stokes wedge defining the original integral 
transform. This is in complete analogy with the following extension of Hermite polynomials: 
He^ a] (y) = a~ n/2 He\^\(p/</a) = e~ ad ^ 12 if n ; meaning that we can choose a — i\x and talk in 
terms of HeJ^ ] (y). 

As mentioned before, we want to represent these results in terms of Fox's H-function. In 
order to do so, we will use the following conversion table: 



He 2n ((/0 = ( 
= ( 

He 2n+1 ((/0 = ( 
= ( 



n! 2 

-1)" ^ 4^ Hj'J [(1 + n, 1); (0, \, 1); </> 2 ] ; 
n! r(— n) • 2 

(2n + l)! o , 

-IT- —^^i-n; 3 -;? 2 ); 

n\ z 

(2n + l)! r(|) 

-1)" H ' [(1 + n, 1); (0, 1); </> 2 ] ; 

n! r(— n) ' 2 



where ^ is Rummer's confluent hypergeometric function, r is the Gamma function, and 
r(i) = V^andr(|) = ^F/2. 
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With this in mind, it is not difficult to write the FPI as an H-transform, 

FPI SDE 

= j>^ 2/ V^ = ae ij2/2 » ; where M = ) ; (20) 

= a \); (0, 1); J] . (21) 

In this particular case, it is possible to reinterpret (20) in terms of an LCT with b = 0, following 
(6). All we need to do is choose 0[J d] =t and, 

where co — i2n\x. This means that i2f M / [J] is a Gaussian Transform of width — i co. 

Now, we can use (17) parametrized by (19) in order to construct Coherent States T£ [J], 
parametrized by \x and (p = ((/?). 

The general expression for T^lJ] is given in terms of the matrix A = (jj) diagonalizing 
M, A" 1 M A = diag(A x , A 2 ), in the following fashion: 

00 

TJJ) = (C A 5J(J) = vTe^ 2 ^ J 5( V - ^e^^e-^^dip ; (22) 

— 00 

where the notation is the same as in (1), (2), (3): that is, the coherent states are the Linear 
Canonical Transforms of the Dirac-delta function parametrized by the diagonalization trans- 
form, i.e., localized states defining each quantum phase (representing each integration cycle 
and Stokes' wedge). 

For this case of ours, the coherent state is given by, 

[J] = y/^ie ij2/2 ^ e in ^ e~ 2ni ^ J ; (23) 



where = M y 4 , and A^ =±2 = =Fl, A^ =±2/ = ^i-iV2. 

The plots below depict T^[J] for /x = ±2 and /i = ±2i: the graphs on the left are the 
wrapping over the Riemann Sphere of the ones on the right. 
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Coherent State plots for different combinations of the parameters /i and (/? - The plots 
on the right have 3t(T^~^Q ,±2l [J]) as the x-axis and 3(Y^~j^ ,±2l [J]) as the y-axis, 
with variations of color (phase 9) and brightness (magnitude p) to indicate z = p e l ° = 
T^j^^fJ]; plots on the left are simply their wrapping over the Riemann Sphere. 



6.2. Cubic Action 



At this point we can try and study a slightly non-trivial model, given by a DO-brane with 
S[if] = (p 3 /3. It's Feynman Path Integral and Schwinger-Dyson equation are shown below: 
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(5/ - J) 3? [J] = => 2f[J] = a Ai[J] + b Bi[J] ; (25) 

where the contour 6 over which the Path Integral is well defined is along the lines of tp such 
that its arguments are either or ±2 tt/3, and Ai and Bi are the Airy functions, that can also 
be expressed in terms of Fox's H-function as 3 2/3 Ai(y/3 1/3 ) = H}J((|, \)\ (0, 1); y/3 1/3 ). 

The Bi function can be expressed in terms of the Ai function via a combination of different 
integration cycles: Bi[J] =e~ in/e Ai[J e~ 2ni/3 ]+e in/6 Ai[J e 2ni/3 ]. That is, we scale and rotate 
the argument. The combination of the two integration cycles can be understood as follows: the 
pre-factor e ±in/e is a critical value for the scaling, and g = e ±2ni/3 is the appropriate Stokes' 
factor associated to crossing each cycle (upon which we cross a cycle and pick up another 
Stokes' factor). 

Now, we want to choose an appropriate M and Q[tp] in order to parametrize our Path 
Integral along the lines of (5). Let us choose M Ai = ~ 2 n ) and 0[(p] = e iLp3/3 . 
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M[J( 
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M[J( 



-2ni/3-\ 



The shaded regions are Stokes' wedges and 
represent zones of convergence: cycles be- 
ginning and ending within these regions 
will yield meaningful integrals. G 1 and C 2 
are particular cycles giving rise to the two 
independent solutions indicated with their 
respective Stokes' factor. 



As we mentioned before, M Ai can be decomposed in the following fashion: M Ai = M DT (— 2 ri) 
M FT M GT (0), i.e., we use a composition of a Gaussian Transform with a Fourier Transform and 
with a Dilation. However, because we now have an Q[ip] = e lip 13 that is non-trivial, it is easier 
to read the analytic structure of our Path Integral straight from its representation in terms of 
its associated H-function. 

As a passing comment, let us just note that (25) can be easily generalized in two different 
ways: 



{df±g 2 J)2? MAi [J} = 0; 
(5 J 3 -4J5 J -2)^ Mi [J] = 0; 
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that are finite at the origin, and where the first generalization can be understood in terms of 
a scaling by a factor of g 2/3 , while the second generalization can be interpreted in terms of 
completing the squares — as expected, these solutions can easily be expressed in terms of H- 
functions (via dilations of its variable or by taking powers of it), and are also parametrizable by 
distinct M's. That is, completing the polynomial potential such that we have an Action along 
the lines of S[ip] = a(p 2 /2 + b(p 3 /3, the coefficient of the highest power ("top coupling") 
determines the analytical structure. As b — > 0, a plays a more and more relevant role: only a 
cycle contained in both regions would yield convergent solutions for both powers — a cycle in 
non-matching convergence regions yields a divergent Sf, while an integration cycle that can 
be deformed from one convergence region to another yields a vanishing 3f (analogous to the 
accumulation of Lee-Yang zeros). 




Therefore, the quantum phases of this system are given by the Ai and Bi functions, and 
their integral representation corresponds to the appropriate Path Integral associated to each 
quantum phase. Moreover, their respective coherent states are given by, 

T^[J] = y^e" j2/2 e-^oV° J ; (26) 

where (p is a constant (something like (p = (Ai|(/?|Ai) or (p = (Bi|<p|Bi), for each quantum 
phase). 

Below you see plots of (26) for (p = ±1: the plots on the right side are such that the 
complex argument (phase) is shown as color (hue) and the magnitude is show as brightness; 
the plots on the left side are nothing but the wrapping of the former over the Riemann Sphere 
(keeping in mind that its group of automorphisms is given by SL(2, C)/{±1}). 
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Coherent State plots for each quantum phase. The plots on the right have 
9t(T (/?o=±1 [J]) as the x-axis and 3(1^ =±1 [J]) as the y-axis, with variations of color 
(phase 9) and brightness (magnitude p) to indicates = pe l9 = =±1 [J]; plots on 
the left are simply their wrapping over the Riemann Sphere. 



6.3. Quartic Action 

The quartic model is that of a DO-brane with S[<p] = g (p 2 /2 + (p 4 /4, where g = [Ji 2 /X just 
for convenience's sake, i.e., g — > +oo corresponds to the weak coupling regime, while g — > 
represents the strong coupling one, and g — > — oo represents the "broken symmetric" one. 
Accordingly, the Path Integral and Schwinger-Dyson equation for this model are given by, 

3T[j] = Je^V24V/4) e iJ</> < oo ^ e : Arg(y) = {0, ±n/2} ; (27) 

e 

(5/ -gdj-J) 3f [J] = => fSf [J] = a t/[g, J] + 6 V[g, J] + c W[g, J] ; (28) 

where the cycle G over which (27) is well defined is such that the arguments of <p are either 
or ±7r/2, and [/, V and W are the Parabolic Cylinder Functions, that can also be cast in terms 
of Fox's H-function as follows: 
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y,ig; v) = PCyl even (g; y) = e^ 2/4 h£ ((^l, 1); (0, 1, 1); ; 

v 4 V 

J 2 (g; V) = PCyl odd (g; y) = y e~^ 4 h£ ((^, 1); (0, - 1, 1); f ) ; 

v 4 y 



where 17", V, W are particular combinations of PCyl even and PCyl odd satisfying the boundary 
conditions of (28) or, equivalently compatible with the integration cycle chosen in (27) (in 
the sense of yielding a well-defined measure assuring the convergence of the integral). Thus, 
the Path Integral of this model is given by the integral representation of Fox's H-function, 
(7), making it clear how the coupling constant g determines the singularity structure and 
integration cycle of the particular Path Integral in question. 
We can now appropriately parametrize (27) using M PC x = ( I_ Q n ) and Q[ip] = e lLp4/4 : 

2 71 

in this way the quantum phases of our model are given by the Parabolic Cylinder Functions 
with respect to the physical parameter g — varying g will change the ground state across 
Stokes' lines due to quantum fluctuations. As before, M PCyl can be understood in terms of a 
composition of a Gaussian Transform of width co = g/2in with a Fourier Transform and with 
a Dilation of A = -2 n: M PC l = M DT (A = -2 ri) M FT M GT (co = g/2n). 




The shaded regions represent convergence 
contours for the following terms, respectively: 
red <^ quartic, blue <^ cubic, and yellow <^ 
quadratic. 



The coherent states for this model can also be computed, yielding, 



T* d [J] = ^exp{, I7 =!=- J } expfi^C^i} exp{-2i Vo ; (29) 

where (p is a constant (a la (p = (U\<p\U), for each quantum phase: U, V,W), and g is the 
model's coupling constant. 
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Below you find plots of (29) for g = +2, 0, —2 and (p = +2 i, 0, +2, respectively. The plots 
on the right side show the phase as color (hue) and the magnitude as brightness. The plots on 
the left side are the wrapping of the former over the Riemann Sphere (keeping in mind that 
its group of automorphisms is given by SL(2, C)/{±1}). 




Coherent state plots for each quantum phase. The plots on the right have 
^(^^"=+2) 0+2^ as x_axis ancl 3(T^~^2°/q^ 2 ) as y-axis, with variations of color (phase 
9) and brightness (magnitude p) to indicates = pe l ° = T^~^ 2 /~^ 2 . Plots on the left 
are simply their wrapping over the Riemann Sphere. Note that the values for g and (p 
are obtained from their relation established via the Schwinger-Dyson equation: <p = 
or if = ±^/—2g , where we have chosen the positive branch of the square root. 
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It is very interesting to note that *+ 2£ — > T 8 ^ = ^ 2 through a rotation of n/2 (J — > i J and 
we also pick up a Stokes' factor, the multiplier e~ l8n ), which can be interpreted in terms of the 
rolling of the Riemann Sphere around the base manifold (which is just a point in this case). 
Moreover, the fixed-point locus of our Action is simply given by the set {0,±^—2g }. Both of 
these together imply that the Gromov-Witten invariant for this model is given by the sum of 
the stationary phase approximation (viz. Atiyah-Bott fixed-point theorem) of the Paht Integral 
computed at each point of the fixed-point set, for different values of the coupling constant g: 
this brings the asymptotic behavior of U, V, W to the forefront, making the Stokes' phenomena 
(wall-crossing) associated to this model quite relevant in regards to the asymptotic behavior 
of each quantum phase. 

The Stokes' Phenomena mentioned above is associated to understanding the differential 
operator defining the Schwinger-Dyson equation (28), (dj 3 - g dj — J), as a connection on 
the Riemann Sphere, forming the basis of an interpretation in terms of the Riemann-Hilbert 
correspondence, which can be thought of as describing the Schwinger-Dyson equation (28) in 
terms of the monodromies of its solutions. The underlying geometry is sometimes called "wild 
geometry" (in analogy with "wild ramification" in arithmetic): the [matrix of the] connection 
cannot be reduced to [a matrix] having logarithmic poles, which is in turn related to the fact 
that (28) is represented via Fox's H-function, and not some form of polylogarithm (which is 
very common in models dealing with MHV amplitudes). 

Both of the observations made in the two paragraphs above are intimately related to the 
topic of nonlinear Fredholm theory [24], in fact, with the nonlinear Fredholm theory of the 
differential operator defining the Schwinger-Dyson equation (28), describing the deforma- 
tions of the given pseudoholomorphic curves (determined by the different solutions to (28)). 
As such, as long as the Path Integral representation of each solution, (27), can be understood 
as a Fredholm map, these observations are also related to the construction of Seiberg-Witten 
invariants. 



6.4. Non-Polynomial Actions 



For the case of non-polynomial actions, we can perform the same construction as for polyno- 
mial ones, just the integration cycle may be a bit more intricate. For example, if S = [3 cos 6 
and 2f[J,J] = J^e^ cos0 e JeI e Je ~ l 7)9, the two relevant integration cycles are G 1 and G 2 > 
where = J e - J e2 , as depicted below, 
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7. Dimensional Extensions 

Firstly, it is important to realize and clarify a few different facts: 

Airy Property: the so-called "Airy property", as defined by [31], needs to be satisfied if we 
want to dimensionally extend the models above; 

Topological Field Theory: our DO-dimensional examples above can be understood as TFTs 
in the sense of Atiyah (the spacetime metric does not appear anywhere in these theories, 
and there is no real dynamics or propagation); 

S-model analogy: when we extend our field if from being scalar-valued to being, e.g., 
Matrix-valued, we enter the realm of S-models, in the sense that the space of Matrices 
is considered the target space of our model — the same argument holds true if rather 
than Matrix-valued we have Tensor-valued, ifie or graded (SUSY) ifie Algebra-valued 
fields, or even more generic target spaces; 

BFSS/IKKT analogy: the DO-branes of our examples can be understood as being in the in- 
finite momentum frame (light-cone frame) and, as such, as a random matrix model 
(extending the fields from scalar- to matrix- valued) where the N — > oo is conjectured to 
be equivalent to M-theory — in this sense, matrices of finite N can be thought of as a 
discretization scheme; 

Causal Dynamical Triangulation: Tensor- valued fields are a generalization of Matrix mod- 
els and correspond to dynamical triangulations, [23], and, as such, the results previously 
obtained are relevant to the field of CDTs; 

Group Field Theory: GFTs were also developed as a generalization of Matrix models, where 
the fields are valued in an appropriate if ie Algebra (possibly graded, SUSY) — once 
again, our results are thus relevant to this field as well; 
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Boundary Conditions: the boundary conditions of the Schwinger-Dyson equations (resp. 
integration cycles defining the Feynman Path Integrals) are required to not be affected 
in the process of taking the N — > oo limit, meaning that the quantum phases of the 
system should remain well-defined within their respective Stokes' wedges (topological 
Dp-brane). However, as N — > oo, we have more and more points satisfying the boundary 
conditions, eventually yielding themselves a Dp -brane in their own right. Thus, we can 
revert our reasoning and say that such Dp -branes built out of the suitable boundary 
conditions are responsible for the quantization of our system, in analogy with [18]. 

Therefore, we have the following general picture. We start with a DO-brane model of fields 
that can be either scalar-valued (e.g., M, C), yielding a topological field theory, or valued in 
more general target spaces (Vector, Matrix, Tensor, ifie or graded ifie algebra, etc), yielding 
a topological S-model. 

In any case, the quantization is done via a generalization of the Feynman Path Integral given 
in terms of a Linear Canonical Transform, which is parametrized generically by a matrix M G 
SL(2,C) (keeping in mind that SL(2, C) is a 3-dim c manifold parametrizing our model). This 
parameter M essentially labels the allowed integration cycles that render the Linear Canonical 
Transform finite and well-defined. Furthermore, it is possible to combine different cycles using 
the rule that M 3 = M 2 'M 1 : if M 1 and M 2 are both labels of allowed integration cycles, so is 
M 3 , i.e., 



= _ ie inJHd 2 b-Hd 1 b-') x (30) 
x e -2i7rJ^(b 2 - 1 + b- 1 ) e i7r^ 2 (a 2 b- 1 +a 1 b- 1 ) e iS[^] rp^ % 

This composition rule is analogous to the gluing of bordisms together, following Atiyah's defi- 
nition of a TFT, in the sense of [ , Definition 1.1.1 and 1.1.5, Example 1.1.9]. This is relevant 
for dimensional extensions for the reason that if we want to combine several DO-branes into 
a Dl -brane, certain relations among the boundary conditions of their respective Schwinger- 
Dyson equations must be satisfied in order for this dimensional extension to be "stable" in 
some suitable sense. These relations are more easily encoded in terms of the algebra of the 
parameters M. 

Thus, if we want to dimensionally extend our O-dimensional systems into bona fide n- 
dimensional ones, we need to first pick an integration cycle, i.e., we need to pick an M labeling 
the particular boundary conditions of the Schwinger-Dyson equations, determining which so- 
lution we are about to extend. Then, we can proceed in two ways, 

BFSS/IKKT inspired: Interpreting the model as a S-model (where scalar-valued systems are 
understood and extended to Vector- valued linear S-models), taking the large-N limit 
and Lorentz boosting it to a desired frame (i.e., away from the infinite momentum 
frame); & 

Many-Body System inspired: Constructing a discretized version of the Feynman Path Inte- 
gral out of several copies of our system, extending the class of allowed paths to include 
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not only Brownian paths, but also Levy flights, and even more generic ones, using the 
full power of the Linear Canonical Transformation. 

These two processes, in principle, should give us the same resulting dimensionally extended 
theory Analogously, we can revert the above arguments saying that given a D-dimensional sys- 
tem, we can define its quantization through the implementation of the appropriate boundary 
conditions of its respective Schwinger-Dyson equation (resp. integration cycles of its Feynman 
Path Integral) in terms of suitable Dp-branes: the solutions to the quantum equations of mo- 
tion are only fully determined when we implement its compatible boundary conditions, which 
is achieved through certain Dp-branes that implement them. 



8. Summary and Outlook 

We have shown how to extend the notion of the Feynman Path Integral in a way to include 
more general paths than the original formulation allows (e.g., Brownian paths, Levy flights, 
etc). This extension comes labelled by a matrix parameter that effectively determines the kind 
of path being used. 

Furthermore, we have established the integro-differential nature of the quantization prob- 
lem, associating the boundary conditions of the Schwinger-Dyson equations to the integration 
cycle of the Feynman Path Integral. In this fashion, the matrix label previously mentioned has 
the job to perform the bookkeeping of all possible allowed solutions, a fact that is reflected by 
the existence of more than one integration cycle guaranteeing the convergence of the Feynman 
Path Integral. 

A relevant property of Linear Canonical Transforms that will be important to us is their eigen- 
functions, forming a complete orthonormal set — in analogy, for instance, with the Fourier 
Transform's eigenfunctions, 

2 l/4 

= ^e~^ 2 He n [2jWip] ; (31) 
Vh! 

where the Hermite polynomials are given by, 

He n [ V ]=(-l)"e^ 2 ^e^ 2 / 2 . (32) 

These will give us information about the coherent states associated to the Linear Canonical 
Transform and each of its components (Fractional Fourier, Hyperbolic, Bargmann, etc). 

As a consequence, it is possible to stratify 1 the space of functions of if (with some suitable 
decaying conditions, determined by the boundary conditions of our problem) according to the 
eigenvalues of its Linear Canonical Transform. The choice of Hermite functions in the Fourier 
Transform case, (31), is convenient because of their localization properties in ip and J, thus 

X A stratified space (or filtration) is one that has been decomposed into subspaces called strata, that are required 
to fit together in a certain way. Stratified spaces provide a setting for the study of singularities via Morse 
Theory on manifolds with boundary and manifolds with corners (e.g., orbifolds). 
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allowing for a more general transform to be introduced and used in field-source analysis (i.e., 
quantization) via Path Integrals or Schwinger-Dyson equations, namely the Linear Canonical 
Transform. 

For example, thinking of the Path Integral as a Sum over Paths (considering the space of 
paths suitably stratified according to the problem at hand), the Path Integral as a Fourier 
Transform implies these paths are Brownian; the Path Integral as a Fractional Fourier Trans- 
form implies these paths are Levy flights; and so on. In this sense, we are trying to formalize 
the notion of a "transform" between two [fundamental] groupoids 2 , the space of paths, if, and 
the space of sources, J, where the Path Integral plays the role of this transform. 

On top of the explicit analysis of three different DO-brane examples, we also showed how to 
dimensionally extend such models in two different ways. This enabled us to reverse the line 
of thought and think of the quantization of a system as being determined by the Dp-branes 
that act as boundaries to it. 

As for the connection with wall-crossing and crystal melting, we clearly stated the role 
that Stokes' Phenomena plays in our models, showing how the asymptotic structure of each 
solution is relevant in determining the allowed integration cycles of the appropriate Feynman 
Path Integral representing the system. Along these lines, we used an intrinsic property of the 
Linear Canonical Transformations, namely its symmetry with respect to SL(2,C), plotting the 
relevant objects over the Riemann Sphere, to make a connection with some dualities (S -duality, 
T -duality) and to establish the modular character of the Feynman Path Integral representation 
we used. 

Finally, we have determined the appropriate representation of the Feynman Path Integral of 
each of our examples in terms of Fox's H-function, which is a Mellin-Barnes transform. Thus, 
we brought to the foreground discussions regarding the singularity structure of our models, 
showing that not all systems can be described in terms of Polylogarithms (as happens in the 
case of MHV amplitudes) : more complicated systems require a more elaborate description, 
where the use of the H-function can be particularly handy. This fact is associated with the 
nonlinear Fredholm theory of the Schwinger-Dyson operator in question, which in turn has 
something to say about the construction of Seiberg-Witten invariants in each of our systems. 

To illustrate this point further, we will use some Cerf theory, which is a way to study families 
of smooth functions on smooth manifolds, their generic singularities, and the topology of the 
subspaces these singularities define, as subspaces of the function space. Essentially, two Morse 
functions can be approximated by one that is Morse at all but finitely many degenerate times. 
The degeneracies involve a birth/death transition of critical points. Cerf theory provides a 
framework for analyzing transitions between Morse functions; it is the study of the positive 
co-dimensional strata of the space of smooth functions. 

For example, / t (x) = x 3 + tx is Morse for t > (no critical points) and for t < (2 
critical points), while t = represents a birth/death transition of critical points, i.e., the 
attachment of a 1-handle connecting the manifold f t>Q to f t<0 . Something similar happens 
with f t {x) — x 4 + which is Morse for t > (1 critical points) and for t < (3 critical 

2 A groupoid generalises the notion of a group in several equivalent ways, e.g., an oriented graph, or a category 
in which every morphism is invertible, etc. The fundamental group measures the 1 -dimensional singularity 
structure of a space; for studying "higher-dimensional singularities", the homotopy groups are used. The 
fundamental groupoid, rather than singling out one point and considering the loops based at that point up to 
homotopy, considers all paths in the space (up to homotopy) . 
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points), while t = indicates the birth/death transition of critical points, i.e., the attachment 
of the handlebody connecting the manifolds f t>0 and f t<0 . 

The plots below highlight this construction for two relevant Potentials, the cubic and quartic 
ones: V}(0) = 4> 3 /3- J 4>, = <£ 4 + g<£ 2 . The birth/death transition, at J = and g = 

(marked as "critical points" in the plots), represents the addition of extra handles. We can 
think of this in terms of constructing the Potential manifold via the attachment of handlebod- 
ies. In physics parlance, we can give the fluxes that determine the handlebody decomposition 
of a manifold (in particular, the vacuum manifold, the moduli space) : once the handlebodies 
determined by the appropriate fluxes are known, the desired manifold is determined — this 
is known as the Heegaard diagram of a manifold. Using the help of the contour (equipoten- 
tial) lines on the plots, we can determine the set of flow lines outgoing or incoming a cer- 
tain critical point (birth/death transition), thus defining ascending and descending manifolds 
(analogous to the ones computed in [36]). Finally, the intersection of these ascending and 
descending manifolds marks the Stokes' lines determining the boundaries between adjacent 
Stokes' wedges. These intersections, called "handle slides" , are the Dp-branes establishing 
the desired bounary conditions for our model. 
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Plots of Vj(4>) = 4> 3 /3 - J 4> and V g {<t>) = 4> 4 + g (p 2 for ranges of the parameters J 
and g, showing the "critical point" where the birth/death transition of critical points 
between two Morse functions happens, i.e., the point where the transition handlebody 
should be attached. The points where J, g = represent the transition between J,g > 
and J,g < 0, where an extra monomial gets added to the top-coupling term. 



At this point we can make some comments about similar approaches and constructions 
available in the literature. 



1. Multi-cut matrix models (which can be extended to Tensor models and thus Group Field 
Theory) have been a research topic for a couple of decades now, [33]. However, their 
associated Stokes' phenomena, and their non-perturbative completion, are much more 
modern features. Their connection with our approach has only become clear in the past 
two years or so. Here is a rough summary of the approach, highlighting the common 
points. 

°s> Riemann-Hilbert and Deligne-Simpson problem: given a system of differential 
equations on the Riemann Sphere, we want to reconstruct the system from its 
monodromy data, [11], i.e., we want to solve the inverse monodromy problem 
in order to obtain the suitable Stokes' multipliers and factors. In our case, the 
system of differential equations is given by the Schwinger-Dyson equations, while 
its monodromy is established by the boundary conditions we choose, i.e., by the 
fall-off properties of the fields in our system. In the case of multi-cut boundary 
conditions, these turn out to be the quantum integrable T-systems; and in the case 
of non-critical string theory, this is described by the position of the cuts. 

°^ Non-perturbative amplitude and corrections: the D-instanton chemical potentials 
are the missing information in [perturbative] string theory, [33]. These are equiv- 
alent to the 9 -vacua in [1] which, in turn, are the Dp-branes establishing the ap- 
propriate boundary conditions rendering our Path Integral well defined. Further, 
these 9 -vacua correspond to the appropriate Stokes' multipliers and factors which, 
as seen above, are related to the Riemann-Hilbert and Deligne-Simpson problems 
via inverse monodromy, thus establishing the Stokes' phenomena of our system. 
This enables us to relate different perturbative vacua and study the string theory 
Landscape from first principles, where the relation between different string theory 
vacua in the Landscape happens via gluing the spectral curves non-perturbatively. 

°^ Spectral Curves and perturbative string theory: the resolvent of the multi-cut ma- 
trix model defines the spectral curve, and thus the perturbative correlators, in the 
following way, 



R(x)= ( — tr 1 \ ; 

\n x-m/ 
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where, 

r 

e - Nuvm dM; 



r 

2? = 

J 



and the perturbative correlators are given by, 

R n (x 1 ,...,x n )=(Y\-tr — 



\ i = l I 

i.e., a N-body problem with potential V (and eigenvalues The resolvent op- 
erator allows us to determine the position of the eigenvalues and the eigenvalue 
density [33]. The non-perturbative corrections are given by the D-instantons, i.e., 
by the 9 -vacua parameters that we compute via the Stokes' factors of the model. 
The geometric meaning of these 9 -vacua is to determine the position of the "eigen- 
value cuts", i.e., the determination of the Stokes' wedges. 

Quantum Integrability: Stokes' phenomena of the Schwinger-Dyson Equations sat- 
isfy the T-systems of quantum integrable models, linking the spectral analysis 
of differential equations and integrable models, establishing a Differential Equa- 
tion/Integrable Model correspondence, [33,34]. 

PT -symmetric Quantum Mechanics started as a proposal to broaden the usual formula- 
tion of Quantum Mechanics by relaxing the requirement of Hermiticity of the Hamilto- 
nian, [35]. Rather than using Hermiticity in order to guarantee that the Hamiltonian 
has a Meal spectrum, it was required that the Hamiltonian satisfied the weaker condition 
of PT symmetry, i.e., that the Hamiltonian be invariant by parity reflection and time re- 
versal. Finally, this generalization is usually studied via a Complex deformation of the 
harmonic oscillator, H = p 2 + x 2 (i x) e , where eGK, exhibiting two phases: when e ^ 
and when — 1 < e < 0; the phase transition occuring at e = 0. 

This scenario can be understood in terms of appropriate integration cycles for the Path In- 
tegral: effectively, the integration cycle is related to e in the following way. Let us choose 
e = 1 to represent the PT -symmetric phase, and e = — 1 for the broken-PT one. When 
e = we have a simple harmonic oscillator (or, in QFT terms, a free theory). Now, we 
can associate a Matrix model to each one of these phases, the PT -symmetric one being 
given by V symm (X) = iX 3 , and the broken-PT one being V bsymm (X) = -iX 2 X _1 . More- 
over, as long as the broken-PT model also has detX ^ 0, we have V bsymm (X) = -iX. 
Thus, from their polynomial form alone we can already see that the same integration 
cycle will not guarantee the convergence of both models. In this sense, PT-symmetry 
can be associate with the suitable cycle guaranteeing the convergence of each respective 
Path Integral. 

Now, following the guidance of our examples, let us define a model without a definite 
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PT-symmetric behavior, given by V(X) = V symm (X) + VbsymmOT = * (uX 3 - v X), where 
we have introduced the modular parameters u, v e C, and where the moduli space is 
given by M = CP 1 , once only the ratio u/v is relevant. This is in complete analogy 
to what was done in subsection 6.2. As such, u (i.e., the top coupling) controls the 
analytic structure of the theory, and as u — > 0, v plays an ever more relevant role, and 
only an integration cycle contained in both Stokes' wedges would yield a convergent 
solution for both powers. As before, an integration cycle in a non-matching region of 
both Stokes' wedges yields a divergent fSf , and an integration cycle that can be deformed 
from one wedge to another yields a vanishing fSf , analogous to Lee-Yang zeros (see the 
plots in subsection 6.2). In fact, we can recast this model as V(X) = igX 3 , where 
g = —u/v G CP 1 , and the linear term, iX, is absorbed in the source term, JX — > JX, 
where J = J + 1. Thus, we end up with a complete parallel to the example of subsection 
6.2, scaled by a [power of the] factor g. Therefore, we can choose a particular linear 
combination of the cycles we obtained in subsection 6.2 in order to endow one of the 
solutions of this theory to be PT-symmetric, while the other solution will have broken- 
PT symmetry. This will be reflected in the allowed ranges of the coupling g, just like the 
Ai and Bi solutions in subsection 6.2 depend on the integration cycles themselves. As 
such, we have a phase transtion at g = 0, and symmetric and broken-symmetric phases 
for 31(g) > and 31(g) < 0. 

3. The Complexification of the Path Integral of Quantum Mechanics, the analytic contin- 
uation of Chern-Simons theory, as well as that of Liouville theory, have recently been 
shown to relate branes in a two-dimensional A-model with possibly new integration cy- 
cles for their respective Path Integrals, [36]. The novelty introduced in these works 
is that a new integration cycle for the Path Integral is possible considering Complex- 
valued paths that are boundary values of pseudoholomorphic maps, giving a middle- 
dimensional cycle in the loop space of a symplectic manifold (classical phase space), 
which is one of the main ideas in Floer cohomology [36]. 

A pseudoholomorphic curve (or J-holomorphic curve, where J is a suitable almost- 
Complex structure, J 2 = 1) is a smooth map from a Riemann surface, S, into an almost- 
Complex manifold, satisfying the Cauchy-Riemann equations, [40]. Moreover, if £ is 
compact (e.g., Riemann Sphere, CP 1 ), there is a nonlinear Fredholm theory describing 
the deformations of the given pseudoholomorphic curve, i.e., these deformations are 
parametrised by a finite-dimensional moduli space which, in turn, will be smoothly de- 
formed by variations in J, [40]. This can be useful in a strategy where J-holomorphic 
curves are used as sigma-models, a : S — > (M, J), since the properties above can be 
used to determine the moduli space of the theory. Further, if J is compatible with a 
symplectic structure on M (think Phase Space), this provides us with an integrability 
condition extending this framework to a well-defined global theory (symplectic topol- 
ogy), [40]. In such a case, I = J^a is a topological invariant of a, and controlled by 
topological data. In physics parlance, if J is compatible with the Poisson bracket of our 
sigma-model's Phase Space (a fact that will depend on the boundary conditions of the 
Schwinger-Dyson equations of the model), the Action of the theory is a [homotopy] in- 
variant. As such, a yields numerical invariants: the Gromov-Witten invariants, opening 
the door to the use of Floer homology, Fukaya category (mirror symmetry, Khovanov 
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homology), and, in four dimensions, to Seiberg-Witten invariants, [40]. 

In [36], the cr-model mentioned above is a topologically twisted A-model, where the 
target space is the Complexification of the original classical Phase Space. The integration 
cycle is described by an A-brane called a coisotropic brane, establishing a relationship 
between the A-model and quantization. In this sense, the Path Integral localizes in the 
new integration cycle (s) of the A-model. 

This is not different from what we presented in this paper (particularly in the case of 
fields valued in a graded Lie algebra, so we can incorporate SUSY to our model and can- 
cel eventual bosonic determinants appearing through the implementation of the bound- 
ary conditions via delta-functionals, on top of being able to use the superpotential in 
order to make things a bit simpler), provided we never cross the Stokes' lines delim- 
iting a certain Stokes' wedge. Otherwise, we have to account for the Stokes' factor 
(hyperasymptotics) associated to this crossing: the bottom-line is the same as in [36], 
the single-valuedness of the integrand, i.e., the choice of the n th root of the canonical 
bundle is an important part of the framework guaranteeing its global existence. Further- 
more, the coherent states, T, that we computed in the examples are Poincare dual to 
their respective defining integration cycles. 

Finally, this is a good point to bring forward the similarities between the present paper, 
[19], and [36] . For example, integral discriminants are an intrinsic property of the Action 
S, as shown below, 



e"V^if^d> ; (33) 



where the simplest cases are just powers of the algebraic discriminant, which determines 
only the singularities of J n \ r . The theory of integral discriminants is related to invariant 
theory, and can be viewed as one of the branches of nonlinear algebra, [19, and refer- 
ences therein] . 

The connections are not difficult to be seen: equation (33) can be understood in terms 
of a Vector field, or we can think of the ip^s as n matrices in a multi-Matrix model, or 
we can interpret the fields in terms of a Tensor model (Group Field Theory). It is very 
interesting to note that [19] expresses the integral discriminants in terms of generalized 
hypergeometric functions, which is in complete analogy with our use of Fox's H-function 
to describe the Path Integral, given that the H-function is a much broader generalization 
of hypergeometric functions, bringing the singularities of the theory to the foreground 
(in parallel to the relation between integral and algebraic discriminants). 

Moreover, [19] uses Ward identities in order to specify appropriate integration cycles 
for the integral discriminants at hand. This is the same role played by the Schwinger- 
Dyson equation (s) in our framework. It is interesting to see the use of a diagrammatic 
technique representing the tensor fields of our model, [19]: this gives us a graphical way 
to compute integral discriminants (Path Integrals), which we have not explored in this 
paper — although it is possible to reconstruct these diagrams from the particular form 
of the H-function representing our theory. 

Further, by /3 -deforming our Matrix model representation (i.e., assuming [3 ^ 1 in the 
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equation below), we can make connections with the Seiberg-Witten prepotentials, 



ar = 



where the Seiberg-Witten prepotential is given by log(2f), [19]. 

It is definitely a very positive reinforcement to see so many equivalences among such 
varied topics and approaches, and the framework we have constructed in this paper. 

4. Complex Action Theory extends the position and momentum operators to be non-Hermitian, 
as well as their eigenvalues to be Complex, allowing a description of a theory with a 
Complex action or of a theory with a Meal action and Complex saddle-points, [37]. The 
basic idea is twofold: 

a) Redefine the inner product of the quantum states in a physically reasonable way, 
such that the Hamiltonian's eigenstates are orthogonal with respect to it, i.e., the 
Hamiltonian has a well-defined Hermiticity under this new inner product — let us 
call this inner product I Q and say the Hamiltonian is Q-Hermitian; & 

b) Suppress the effect of the anti-Q-Hermitian part of the Hamiltonian after long times, 
thus the effect of the imaginary part of the Q-Hamiltonian (the anti-Q-Hermitian 
part) is smoothed out except for an irrelevant constant. 

Compared to our approach, this is the same as realizing that the integration cycle (i.e., 
the particular Stokes' wedge where the Path Integral is defined) determines the allowed 
values of the parameters of our theory (the coupling constants), and that gives us control 
over the asymptotic behavior of our theory: that is, ultimately the selected Stokes' wedge 
determines the parameters of the theory which, in turn, control the fall-off behavior of 
the fields in the theory and thus its asymptotic properties. 

In some sense, this is not very different from the PT-symmetric Quantum Mechanical 
case mentioned above, where the inner product is modified via a charge conjugation 
symmetry, C, such that CPT is ultimately conserved in the theory. In this fashion, as soon 
as we choose a suitable integration cycle retaining PT symmetry, it suffices to choose a 
compatible I c inner product. Thus, this charge conjugation symmetry, C, plays the role 
of the deformation Q mentioned above, [37], defining a C-Hermitian property. This can 
be understood in terms of the extensions of operators, as done in Appendix A, where the 
extension is given by C. In short, the properties of the particular chosen integration cycle 
select the parameters of the model which determine not only its asymptotic behavior 
but also the C-extension, clearly implying that different Stokes' wedges (i.e., integration 
cycles) yield distinct sets of coherent states labelled t c = J e |$)($| 2)$, where 6 is the 
integration cycle associated to the charge conjugation symmetry, C. 

Complex Action Theory is then an extension of PT-symmetric Quantum Mechanics in 
the sense that one is free to pick any of the allowed integration cycles and define the 
charge conjugation symmetry, C, through it. Then, it suffices to define a compatible PT 
symmetry, and build the theory from there. This enables the use of all Stokes' wedges, 
as opposed to the situation in PT-symmetric Quantum Mechanics. 
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5. A new framework for Complex- valued classical fields in the case of quantum field the- 
ories describing neutral particles has recently been proposed, generalizing Osterwalder 
and Schrader's construction of Euclidean fields via reflection positivity, studying scenar- 
ios that ensure and preserve it, [38]. 

The relevant part of this framework is that the neutral fields can be either Meal or 
Complex — the usual distinction between neutral and charged fields does not corre- 
spond with the distinction between Meal- and Complex- valued fields: Complex neutral 
fields allow for some novelty. As made quite clear in [38], the two-point function of the 
fields is the integral kernel of an operator, D, which need not be Hermitian. However, 
two requirements are made: 

a) The Hermitian part of D should have a strictly positive spectrum; & 

b) The transformation D should be reflection positive. 

Charged fields are also treated in an innovative fashion: distinct charged fields are in- 
troduced, not related by Complex conjugation, where charge conjugation acts unitarily 

Reflection Positivity not only gives meaning to inverse Wick rotations, but is also used 
in the analytic continuation of group representations, in understanding the properties of 
the spectrum of the Transfer matrix, and in the theory of phase transitions, [38]. 

The connection to our approach can be established when we are able to use the func- 
tional integral description of the above, which describes the decaying properties of the 
fields of the theory, i.e., their fall-off properties at the boundary. This is equivalent to set- 
ting boundary conditions to the Schwinger-Dyson equations appropriate to our model. 
In this sense, the problem of quantization is reduced to that of establishing suitable 
boundary conditions determining the desired decaying properties for the field content 
of our theory This procedure can be accomplished in different ways, one of them being 
the construction of an appropriate Fourier Transform respecting the constraints of the 
problem, that is then understood as the partition function for the model (in an analogy 
with the Paley-Wiener theorem, [16]). 

Essentially, this can be understood as an extension of PT-symmetric Euclidean QFT, 
where we can choose different combinations of C, P, and T to be held by the Hamil- 
tonian, and the remaining symmetry "corrects" the inner product. For example, we can 
choose a CP-symmetric Euclidean QFT whose inner product is corrected by T, etc. As 
such, this can be understood in terms of J-extensions of J-symmetric operators, [7,8,10]. 

In this sense, this discussion can be reduced to that of PT-symmetric Quantum Mechan- 
ics and of Complex Action Theory above. 

As a final remark, it is worth saying that we tried to understand the process of quantization 
via the singularity structure and poles of the path integral, 2ff[J] 9 appropriately decorated 
by the parameters of the theory, t, and suitably compatible boundary conditions, S. In the 
language of [41], we can loosely say that we tried to study the functoriality of quantization 
via the singularity structure of the path integral: facts on a topological setting (e.g., fall- 
off properties of the fields, asymptotic structure of the path integral) can be translated to 
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an algebraic setting (algebra of observables, classification of the vacua of the theory). The 
representation of the path integral as a Fox H-function encodes the analytic structure of the 
theory in a Mellin-Barnes transform, which can be understood as an appropriate L -function 
(viz. functoriality via poles of L -functions), [41]. 

In particular, as suggested in [41], one can think of Quantization c Deformation: study 
the symmetries of a certain equation of motion to be quantized and deform the Lie algebra 
associated to each symmetry. These should be able to be expressed in terms of their own path 
integral. Deformation and Functoriality then can be seen as a reflection of the wave-particle 
duality, respectively: functoriality ~ particle ~ Heisenberg picture, deformation ~ wave ~ 
Schrodinger picture. 
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A. Extensions of Operators 

Another approach to the material developed above would be in terms of extensions of symmet- 
ric operators, [7]. Denoting the Schwinger-Dyson operator by ^ SD , it is said to be symmetric 



Its self-adjoint extensions can be understood in terms of its boundary values (resp. Dp- 
branes), i.e., the monodromy of the boundary conditions (resp. brane monodromy) param- 
eterize the possible extensions. Analogously, we can study ^ SD via its Cayley transform, 
U G = (^ SD — il)(^ SD + il)~\ where 1 is the unit operator: if U G is unitary then ^ SD is 
self-adjoint. This can be measured through the deficiency indices of ^ SD , which are defined 
as the dimension of the orthogonal complements of its domain and range: if both deficiency 
indices are identical then U is unitary. 

The interesting thing to notice is that different self-adjoint extensions may have different 
spectrum. That is, by varying the monodromy of the boundary conditions (i.e., varying the 
brane monodromy) we obtain different operators 6'^ which may define different theories in 
terms of their particle spectrum. 

This can be further understood in terms of an extension of the resolution of the identity, i.e., 
coherent state quantization a la [11], 



if 



(^SD<£ly) = (<£I^SDV>. 
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where E m= i^ n labels the integration cycle associated to each monodromy generator, and 
e lc "*=w represents the contribution from each cycle. 

This can be further generalized by relaxing the condition that the extension be self-adjoint, 
[10]. Instead, we can talk in terms of J-extensions, where J is a conjugation: (J 4>,J if) = 
{4>, ip). That is, a J-extension of SD is such that J SD J = 0^ <^> (^ SD 4>, J if) = {<p,J SD if). 
For example, <B T-symmetric quantum mechanics falls in this case if we choose J = ( P ( T, where 
we have to modify the inner product slightly in order for J to remain a conjugation (so we can 
still retain CPT invariance, [13]). 

This should not necessarily come as a surprise because of the following observation. We 
can understand the boundary conditions for SD via explicit boundary terms in the Action. 
These explicit boundary terms are implemented via new fields that are subjected to the actual 
boundary conditions we desire. The new fields are related to the original ones: they are merely 
convenient objects used to express the constraints necessary for the boundary conditions of 
SD to be [explicitly] reflected in the Action. Finally, these boundary terms in the Action 
are equivalent to Dp-branes described by these new fields: this is the rationale behind brane 
quantization, [18]. 

Following this construction, the notion of J-extensions described above is translated in terms 
of deformations of the integration cycle for the path integral, done in such a way as to imple- 
ment the particulars of the conjugation at hand. 
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